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_j 1 Abstract 

In this paper, we propose a new language, called AR (Action Rules), and describe how 
various propagators for finite-domain constraints can be implemented in it. An action 
rule specifies a pattern for agents, an action that the agents can carry out, and an event 
pattern for events that can activate the agents. AR combines the goal-oriented execution 
model of logic programming with the event-driven execution model. This hybrid execution 
model facilitates programming constraint propagators. A propagator for a constraint is an 
agent that maintains the consistency of the constraint and is activated by the updates of 
the domain variables in the constraint. AR has a much stronger descriptive power than 
indexicals, the language widely used in the current finite-domain constraint systems, and 
is flexible for implementing not only interval-consistency but also arc-consistency algo- 
rithms. As examples, we present a weak arc-consistency propagator for the all_distinct 
f*"*) , constraint and a hybrid algorithm for n-ary linear equality constraints. B-Prolog has been 

extended to accommodate action rules. Benchmarking shows that B-Prolog as a CLP(FD) 
, system significantly outperforms other CLP(FD) systems. 
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1 Introduction 

CLP(FD), the constraint logic programming language over finite domains, has 
been proved effective for solving a large number of real-life optimization prob- 
lems (Din cbas et al. 19901 IJaffar and Maher 1 994). The key operation employed in 
CLP(FD) is called constraint propagation IjKumar 19921 |Tsang 1993| ), which uses 
constraints actively to prune search spaces as follows: whenever a variable changes, 
i.e., the variable has been instantiated or its domain has been updated, the do- 
mains of all the remaining variables are filtered to contain only those values that 
are consistent with this variable. There may exist different propagation rules for a 
constraint depending on the level of consistency to be achieved. Constraint prop- 
agation has been employed to solve not only constraints over finite domains but 
also constraints over trees, lists, finite sets, floating-point numbers, and many other 
domains ( .Taff ar and Maher 1 994). 
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In early CLP(FD) systems, such as the CHIP system (Dincba s et al. 1988J) . con- 
straints are interpreted rather than compiled. Constraints are first transformed into 
canonical-form terms and are then executed by an interpreter that performs, among 
other things, constraint propagation. The propagation procedure adopted is general 
enough for handling all types of constraints. Learning from the success of compil- 
ing Prolog programs into the Warren Abstract Machine (WAM) (Warren 1983), 
a former research group at ECRC extended the WAM for compiling CLP(FD) 
( |Aggoun and Beldiceanu 1991| ). The CHIP compiler compiles constraints into low- 
level instructions such that different specialized propagation procedures are used 
for different types of constraints. This black-box approach has proved problematic 
because it is too complicated and lacks flexibility and extendibility. The extended 
WAM in the CHIP system (A ggoun and Beldiceanu 1991| ) has over 100 instructions 
for compiling finite-domain constraints alone! 

A language construct, called indexicals, has been quite popular as an intermediate 
language for compiling finite-domain constraints. The language was first proposed in 
dvan Hentenryck et al. 1992) ) and then popularized by ( |Codognet and Diaz 1 996). 
This language is also adopted by other systems l|Carlsson et al. 19"97llSidebottom and Havens 1 996 ) 
An indexical is a primitive constraint in the form of X in r, where X is a domain 
variable and r is a range expression for X. For each indexical, a propagation pro- 
cedure specific to it is used. Indexicals are claimed to be a glass-box approach to 
compiling constraints in contrast with the black-box CHIP compiler. Nevertheless, 
as the delaying mechanism is embedded in range expressions, indexicals are not as 
open as claimed. Indexicals can be used to compile arithmetic constraints, but are 
too weak to be used to program many other kinds of propagators. 

CHR (Constraint Handling Rules) (Fruhwirth 1998) may currently be the most 
powerful implementation language for constraints. It can be used to program not 
only constraint propagators but also constraint reasoning rules. CHR has been 
implemented and integrated with ECLiPSe, SICStus, HAL, and Oz. CHR resembles 
a production system. In CHR, the left-hand side of a rule specifies a pattern of 
constraints in the constraint store and the right-hand side specifies new constraints 
to replace those on the left-hand side or to be added into the store. The left-hand 
side of a rule may have multiple constraint patterns. This feature is helpful for 
reasoning about the constraint store. For example, A>Bk,B>C^A> 
C + 1 is a CHR rule that generates the constraint A > C + 1, which is helpful 
albeit redundant. The strong descriptive power, however, is not offered without 
cost. For CHR, a sophisticated matching algorithm is needed to match constraint 
patterns against the constraint store. Now, constraint solvers implemented in CHR 
are still an order of magnitude slower than constraint interpreters implemented in 
C IHolzhaur and Fruhwirth 1 9991 IHolzhaur et al. 200411. 

This paper proposes a new language, called AR (Action Rules), which can be used 
to program event handling in general and constraint propagation in particular. An 
action rule specifies a pattern for agents, an action that the agents can carry out, 
and an event pattern for events that can activate the agents. An agent behaves in an 
event-driven manner. An agent can be suspended when certain conditions on it are 
satisfied and can be activated when certain events are posted. AR is an extension of 
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delay constructs such as delay clauses HMeier 1994(1 that allows for the descriptions 
of not only delay conditions but also activating events and actions IjZhou 1 998). 
The syntax, operational semantics, and implementation of AR will be described in 
Section 3. 

The focus of this paper is on how to implement various propagators for finite- 
domain constraints in AR. A propagator for a constraint is an agent that maintains 
the consistency of the constraint and is activated by the updates of the domain vari- 
ables in the constraint. In Section 4, we present propagators for binary, non-binary, 
and the global constraint all_distinct. AR is more expressive than indexicals. 
Some of the propagators presented, such as the one for maintaining arc consistency 
for binary equality constraints and the one for maintaining weak arc consistency 
for all_distinct, cannot be implemented in indexicals as efficiently. 

B-Prolog has been extended to accommodate AR and several constraint solvers 
including the ones over finite domains, Boolean, trees, and hnite sets have been 
developed in AR (Zho u 2002(1 . Section 5 compares the performance of the finite- 
domain solver of B-Prolog with GNU-Prolog (GP), a state-of-the-art implemen- 
tation of CLP(FD) ( |Diaz and Godognet 200T| |, and two other CLP(FD) systems: 
ECLiPSe and SICStus. Benchmarking results show that B-Prolog is significantly 
faster than GP and 4-6 times as fast as ECLiPSe and SICStus. 

Readers are assumed to be familiar with logic programming and constraint sat- 
isfaction, but no knowledge about the compilation is assumed. In Section 2, we 
define some preliminary terms and concepts about CLP(FD) and constraint prop- 
agation. Readers are referred to ( |Marriott and Stuckey 1998| l, ( |Hentenryck 1989| ) 
and l|Kumar 1992(1 for the details. A brief description of the implementation of AR 
is given in Section 3, and a more detailed description can be found in ((Zhou 2003 ). 



2 Preliminaries 
2.1 CLP(FD) 

CLP(FD) Pentenryck T989| ) is an extension of Prolog that supports built-ins for 
specifying domain variables, constraints, and strategies for instantiating variables. 
The domains of variables are declared as follows: 

Vars in D 

where Vars is a variable or a list of variables, and D is a list of ground terms 
or a range between two integers l..u. A domain variable is normally represented 
as a Prolog variable with attributes. A CLP(FD) system provides primitives for 
accessing and updating attribute values. 

A CLP(FD) system provides equality (=), disequality (7^), and inequality con- 
straints. In addition, a CLP(FD) system also provides some other constraints such 
as global constraints. The global constraint all_distinct (L) ensures that the el- 
ements in the list L must be all pairwise different. 
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2.2 Constraint propagation 

Constraint Propagation IjKumar 19921 |Tsang 1993| ) is a key operation employed in 
CLP(FD) systems for maintaining the consistency of constraints. The basic idea of 
constraint propagation is to activate the propagators of constraints whenever the 
domains of the variables in the constraints are updated. Propagating the updates 
to other variables may result in the shrinking of the domains of the variables or the 
instantiation of the variables. 

There are different levels of consistency for constraints such as node, interval, 
bounds, arc, and path consistency ( |Tsang 1993| |Marriott and Stuckey 1998] ) . We 
define below three levels of consistency needed in this paper, namely node, interval 
and arc consistency, and define the propagators that maintain them. 

A unary constraint p(X), where X has the domain D, is said to be node- consistent 
if, for any element x in D, p(x) is satisfied. 

y xeD p(x) 

For example, for the equality constraint X = Y + l, when X is instantiated to 3, the 
constraint becomes unary and Y must be instantiated to 2 to make the constraint 
node-consistent. As another example, for the disequality constraint X ^ Y, when 
X is instantiated to 3, 3 must be excluded from the domain of Y to make the 
constraint node-consistent. The propagation rule that maintains node consistency 
is called forward checking. A propagator for a constraint that performs forward 
checking is activated whenever the constraint becomes unary. 

Let C be a linear equality constraint c + ttixl 1 +a 2 x X% + . . . + a n x X n = 
where a, ^ and Xi is defined over the domain Di (i = 1, . . . , n). Let 

9i {Xi , . . . , Xi- 1 , Xi + i , . . . , X n ) = 

-c - oi x Xi - ... - a l -i x Aj_i - a i+ i x X i+ i - . . . - a n x X n 

at 

and I and u be the functions defined as follows: 

l(9i{Xi, ■ • ■ , Xi-i, Xi+i, . . . , X n j) = 

minig^xi, . . .,Xi-i,x i+ i, . . .,x n )\x k e D k , l<k<n,k^i} 

u (gi(Xi, . . . , Xi + \, . . . , x n )) = 

max{gi(xi, . . . ,Xi-i,x i+ i, . . . ,x n )\x k e D k , 1 < k < n, k ^ i} 
The constraint C is said to be interval consistent w.r.t. Xi if: 

VxeDi(l(gi{Xi, . . . , Xi-i, X i+ i, . . . ,X n )) < x < u(gi(Xi, . . . ,X t -i,X i+ i, . . . ,X n )) 

To make the constraint interval consistent w.r.t. Xi, we have to exclude all the 
elements from Di that are not in the range. The constraint C is said to be interval 
consistent if C is interval consistent w.r.t. all the variables. For example, the con- 
straint X = Y + 1, where X and Y have the domain 1..5, is not interval-consistent. 
To make the constraint interval consistent, we have to exclude 1 from the domain 
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of X and 5 from the domain of Y. Propagators for maintaining interval consistency 
are activated whenever a bound of a variable is updated or whenever a variable is 
instantiated. The definition can be easily extended to an inequality constraint. 

Consider a binary constraint p{X, Y) where X and Y arc defined over the domains 
D x and D y , respectively. The constraint is said to be arc- consistent w.r.t. X if for 
any element in D x there exists a supporting element in D y such that the constraint 
is satisfied: 

Similarly, the constraint is arc-consistent w.r.t. Y if for any element in D y there 
exists a supporting element in D x such that the constraint is satisfied: 

The constraint is arc- consistent if it is arc-consistent w.r.t. both X and Y. For 
example, the equality constraint X = Y + 1 (X e {2, 4, 5}, Y e {1--4}) is not 
arc-consistent since there is no element in the domain of X that supports 2 in 
the domain of Y. To make the constraint arc-consistent, we must exclude 2 from 
the domain of Y. Propagators for maintaining arc consistency are triggered when- 
ever changes occur to the domain of a variable. Maintaining arc consistency for 
a non-binary constraint requires examining Cartesian products of the domains 
IjDechter 2003(1 and is thus very costly. For this reason, some CLP(FD) systems 
maintain arc consistency only for binary constraints and many others do not con- 
sider arc consistency at all. 

2.3 Domain variables 

A domain variable is a suspension variable to which there are suspended propagators 
and some other information attached. A domain variable is represented in B-Prolog 
as a record that has the following fields: 



ref reference to the value 

type type of the domain 

min minimum element in the domain 

max maximum element in the domain 

size number of elements that remain in the domain 

ins-cs list of propagators to be executed when the variable is instantiated 

bound jcs list of propagators to be executed when a bound is updated 

domjes list of propagators to be executed when an inner element is excluded 

elms pointer to the bit vector representation of the elements 



where ref refers to the variable itself if the variable is not instantiated, and elms is 
a pointer to a bit vector that represents the elements. When a domain is an interval 
without holes, no bit vector is necessary and elms is a null pointer. 

The following built-in predicates and functions are available on domain variables. 

• dvar(X): Succeeds if X is a domain variable. 



G 



N.F. Zhou 



• min(A) , max(A) : Functions that return, respectively, the minimum and max- 
imum elements of the domain of A. 

• size(X, Size): The size of the domain of X is Size. 

• exclude (X, E) : Excludes the element E from the domain of X. 

• A in D: The new domain for X is the intersection of its existing domain and 
D where D is a set of ground terms or a range l..u of integers. 

A failure occurs when the domain of a variable becomes empty. When the domain 
of a variable becomes a singleton, the variable is instantiated to the element auto- 
matically. 

An event is posted whenever the domain of a variable is updated. For a domain 
variable A, instantiating A posts the event ins (JO, 1 updating the lower or upper 
bound of the domain posts the event bound (X) , and excluding an inner element E 
from the domain posts the event dom( X, E) . Notice that the event bound (A) is not 
posted when A is instantiated and the event dom(A, E) is not posted when cither 
bound of the domain of A is updated. This implies that a propagator that maintains 
arc consistency has to handle not only dom(A, E) events but also bound (A) and 
ins (A) events. 

Each event on a domain variable activates its corresponding list of propagators. 
The event ins (A) activates the propagator list ins-cs of A, bound (A) activates 
the list bound-cs, and dom(A, E) activates the list dorri-cs. 

3 The AR Language 

AR is designed for programming interactive agents. In this section, we describe the 
syntax, operational semantics, and implementation of AR. 

3. 1 Syntax 

An action rule takes the following form: 

Agent, Condition, {Event} => Action 

where Agent is an atomic formula that represents a pattern for agents, Condition 
is a conjunction of conditions on the agents, Event is a non-empty disjunction 
of patterns for events that can activate the agents, and Action is a sequence of 
subgoals. 2 Condition and the following comma can be omitted if no condition is 
needed on Agent. Action cannot be empty The subgoal true represents an empty 
action that always succeeds. An action rule degenerates into a commitment rule if 
Event together with the enclosing braces are missing. Conditions, event patterns, 
and actions are all atomic formulas where the delimiter ',' is used to separate the 
constituents. 

An AR predicate consists of a sequence of rules defining agents of the same 

1 A variable is said to be instantiated if it is bound to another term, possibly another variable. 

2 Subgoals in Action can be any subgoals including those defined by Prolog clauses. 
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predicate symbol. In a program, AR predicates can be intermingled with Prolog 
predicates. In this paper, the term agents is used to refer to subgoals that can 
be suspended and activated, and the term predicate is used refer to both AR and 
Prolog predicates unless explicitly specified. 

All conditions must be in-line tests. 3 In the implementation of AR in B-Prolog, 
the following types of conditions are allowed: 

• Type and mode checking: Predicates like integer (A) , var (X) , and nonvar (A) . 

• Matching: A matching call takes the form of X = Y where one of the ar- 
guments must be a non-variable term at compile time and the other must 
be a variable (again at compile time) that occurs before in the rule. The 
non-variable term serves as a pattern and the variable refers to a term to be 
matched against the pattern. This call succeeds if the pattern and the term 
become identical after a substitution is applied to the pattern. For instance, 
the condition f(X) = Y succeeds if Y is a structure whose functor is f/1. 

• Term inspection: Several built-in predicates including arg/3, functor/3, ==/2, 
\==/2, and n_vars_gt/2 can be used in the condition of a rule to inspect the 
arguments of an agent. The call n_vars_gt (Term, N) succeeds if the number 
of variables occurring in Term is greater than N. 

• Arithmetic comparison: Checks the arithmetic equality (=:=), discquality 
(=\=), or inequality (>, >=, <, and =<) of two terms which must be ground 
at runtime. 

A set of built-in events is provided. 4 As far as programming constraint propaga- 
tors is concerned, an event pattern can be one of the following: 

• generated: The action of the rule is executed when the agent is suspended 
for the first time. 

• ins(X): The agent is activated when an event ins (A) is posted. 

• bound (X): The agent is activated when an event bound (A) is posted. 

• dom(A) and dom(A, E): The agent is activated when an event dom(A, E') is 
posted. Before the action is executed, E is made to reference the element E' . 

A user program can create and post its own events and define agents to handle 
them. A user-defined event takes the form of event (A, T) where X is a suspension 
variable that connects the event with its handling agents, and T is a Prolog term 
that contains the information to be transmitted to the agents. If the event poster 
does not have any information to be transmitted to the agents, then the second 
argument T can be omitted. The built-in post(_B) posts the event E. 

In an action rule, the event pattern dom(A, T) or event (A, T) is not allowed to 
coexist with any other event patterns and T must be a first-occurring variable so 
that when the action of the rule is executed T always refers to the second argument 
of the event. 

3 An in-line call is compiled into instructions that do not invoke any predicates. A test does not 
change the instantiation status of the variables in its arguments. 

4 In the implementation in B-Prolog, built-in events are provided for programming constraint 
propagators, graphical user interfaces, and interactive agents. 
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3.2 Examples 

The following defines an agent that echoes the messages sent to it by event posters. 

echo_agent (X) , {event (X, Message)} => write (Message) . 

The following query, 

echo_agent (Ping) , echo_agent (Pong) , 

post (event (Ping, ping)) , post (event (Pong, pong) ) 

generates two echo agents echo_agent (Ping) and echo_agent (Pong) , and activates 
them by posting two events. The event event (Ping, ping) activates the agent 
echo_agent (Ping) , and the event event (Pong, pong) activates echo_agent (Pong) . 
The following defines the freeze predicate in Prolog-II IjColmerauer 1984jl . 

freeze(X.G), var(X), {ins(X)} => true. 
freeze(X,G) => call(G). 

The primitive freeze (X,G) is logically equivalent to call(G) but the execution of 
G is delayed until X is instantiated to a non- variable term. The agent freeze (X , G) is 
suspended waiting for an event ins (X) when X is a variable. When an event ins (X) 
is posted, the condition var (X) is tested again. If it succeeds, then the action true 
is executed and the agent becomes suspended again. As long as X is a free variable, 
the agent freeze (X,G) is suspended. Only when X becomes a non- variable term, 
can the second rule be applied. 

Consider, as another example, how to implement the following indexical: 

X in min(Y)+min(Z) . .max(Y)+max(Z) . 

which ensures that the constraint X = Y+Z is interval-consistent w.r.t. X. 

'V in V+V (X,Y,Z) .{generated, ins(Y) ,bound(Y) ,ins(Z) ,bound(Z)> => 
reduce_domain(X,Y,Z) . 

reduce_domain(X,Y,Z) => 

L is min(Y)+min(Z) , U is max(Y)+max(Z) , 
X in L. .U. 

The propagator is activated whenever a bound of Y or Z is updated or either one is 
instantiated. The action reduce_domain(X,Y,Z) enforces that the domain of X be 
in the range min(Y)+min(Z) . .max(Y)+max(Z) . The action is also executed before 
the propagator is first suspended so that no preprocessing is needed to enforce 
interval consistency. 



3.3 Operational Semantics 



The operational semantics of AR can be presented as a state-transition system 
as shown in Figure ^ An agent may be in one of the following states: start, sleep, 




Fig. 1. Diagram of state transition of agents. 

woken, and end. When an agent is generated, it enters the start state and is executed 
immediately. A typical agent transits to the end state through the sleep and woken 
states. An agent is said to be floundering if it stays in the sleep state forever. It is 
the programmer's responsibility to prevent agents from floundering. 

When an agent is generated, the system searches in its predicate in textual order 
for a rule whose agent-pattern matches the agent and whose condition is satisfied. 
This kind of rule is said to be applicable to the agent. Formally, an action rule 
"H, C, {E} => £>" or a commitment rule "H, C => B" is applicable to an agent a 
if there exists a unifier 9 such that H6 = a and CO is satisfied. 5 If no rule is found 
applicable, the agent fails. If a commitment rule is found, the agent is substituted 
for the body, 6 and its state is changed to end. 

If an action rule U H, C, {E} => £>" is found for the agent, the agent is suspended, 
transiting from start to sleep. It will stay in the sleep state until it is activated by 
one of the events in E. 

When an event is posted, all the sleeping agents waiting for the event in the 
system are woken up and the event is erased after that so that no agents generated 
later will be responsive to this event. The woken agents are added to the queue of 
active subgoals in some order. It is up to the implementer of the language to use 
a strategy to schedule activated agents. Whatever scheduling strategy is adopted, 
the programmer should not rely on the strategy to guarantee the correctness of 
programs. 

Suppose an agent was put into sleep by the action rule "H, C, {E} => B" and 
was woken up by one of the events in E. After this agent is picked by the scheduler, 
the system tests the condition C again. If it is satisfied, the action B is executed. 
If the action succeeds, the agent is re-suspended. If the action fails, the agent fails 
as well. If the condition C of the action rule does not hold, the system searches for 
an alternative applicable rule for the agent just as for a newly generated agent. 

There is no primitive for killing agents explicitly. An agent never disappears as 
long as action rules are applied to it successfully. An agent transits to the end state 
only when a commitment rule is applied to it. 

5 Notice that since one-directional matching rather than full unification is used to search for an 
applicable rule and in the condition no variable in a can be instantiated, the agent remains the 
same after an applicable rule is found. 

6 A commitment rule is similar to a guarded clause in concurrent logic languages jShapiro 1989} , 
but an agent can never be blocked while it is being matched against an agent pattern. 
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3-4 The Implementation 

The abstract machine of B-Prolog, called ATOAM IjZhou 1996bjl . is extended to 
support agents QZhou 2003(1 . This subsection briefly describes this implementation. 
A more detailed description can be found in IjZhou 2 003 ) . 

3.4-1 The frame structure for agents 

The ATOAM is a variant of the Warren Abstract Machine (WAM) \ Warren 1983|l . 
Unlike in the WAM where arguments are passed through argument registers, argu- 
ments in the ATOAM are passed through stack frames and only one frame is used 
for each subgoal. Each time a predicate is invoked by a subgoal, a frame is placed 
on top of the control stack unless the frame currently at the top can be reused. 
Frames for different types of predicates have different structures. 

Agents are stored as frames on the control stack. The frame for an agent has the 
following slots in addition to those included in a normal frame: 7 

STATE: State of the agent 

EVENT: Activating event 

REEP: Re-entrance program pointer 

PREV: Previous agent in the chain 

The STATE slot has one of those states shown in Figure[T]as its value. The EVENT 
slot stores the last event that activated the agent. The action rule of the agent can 
have access to this event. The REEP slot stores the program pointer to continue 
when the agent is activated. The PREV slot stores the pointer to the previous 
agent's frame in the chain of agents. 

The frames on the control stack comprise three chains, namely the chain of active 
subgoals 8 that are being executed, the chain of choice points, i.e., subgoals that have 
alternative clauses to be tried when execution backtracks to them, and the chain 
of agents in either sleep or woken states. 

Storing agents on the stack facilitates context switching for agents (Zhou 1996a) 
but complicates memory management. With frames of agents on the stack, the 
chronological order of frames is no longer preserved, and therefore a garbage collec- 
tor is needed to collect useless frames on the control stack and run-time checking 
is needed to determine whether the current frame can be reused. 

3-4-2 When and how to invoke agents? 

For the sake of efficiency, events are not checked after each instruction but checked 
at the entry and exit points of each predicate. If it is found that the list of events 

7 A normal frame has the following slots: arguments, FP (parent frame), CPS (continuation 
program pointer on success), TOP (top of the control stack), BTM (bottom of the frame), 
and local variables. A choice point frame has the following additional slots: CPF (continuation 
program pointer on failure), B (parent choice point), H (top of the heap), and T (top of the 
trail stack). 

8 i.e., the frames connected by FP. 
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is not empty, those agents that are waiting for the events are added into the active 
chain and the current predicate is interrupted. After the activated agents com- 
plete their execution successfully, the interrupted predicate resumes its execution. 
The programmer has no control over the order in which agents are added. In our 
implementation, the first- generated- first- served strategy is used. 

At a point during execution, there may be multiple events posted that are all 
expected by an agent. If this is the case, then the agent must be activated once for 
each of the events. If an agent is found to be active already when the system tries 
to add it into the active chain, the sytem makes a copy of it and adds the copy into 
the chain. 

The actions of constraint propagators are to reduce the domains of variables. This 
characteristic is exploited to improve the performance of constraint propagators. 
Some events that cannot lead to the shrinking of any domains are ignored. For 
example, if multiple events of bound (X) are posted at the same time, then only one 
of them needs to be handled, and if bound (X) and ins(X) arc posted at the same 
time, then the bound (X) event is ignored. In this way, many redundant activations 
of rules that do not contribute to the reduction of any domains can be suppressed. 
This optimization is applied to constraint propagators only and not general agents. 

In a CLP(FD) program, constraint propagation is normally intertwined with 
non-deterministic subgoals such as labeling that assign values to variables. If a 
non-deterministic predicate is interrupted by events, then no choice point can be 
created for the interrupted predicate until the activated agents are all executed. 
Consider the following example: 

?-p(X) , X=f (Y) ,q(X) ,write(X) . 

p(X) ,var(X) ,{ins(X)} => true. 
p(X) => X=f (a). 

q(X) :-fail. 
q(X). 

First the agent p(X) is generated, waiting for X to be instantiated. The subgoal 
X=f (Y) posts an event ins(X) after X bound to f (Y). At the entry of q(X), the 
event is detected and p(X) is activated. If there were a choice point created for q(X) 
before p(X) is activated, the binding Y=a given by the second rule of p(X) would 
be lost when fail in q/1 is executed since Y is older than the choice point, and the 
output from write (X) would be f (Y) not f (a). 

4 Programming Constraint Propagators in AR 

The high descriptive power of AR opens new ways to implementing constraint 
propagators. In this section, we implement propagators that maintain node, in- 
terval, and arc consistency for binary constraints, a hybrid algorithm for non- 
binary constraints, and a weak arc-consistency propagator for the global constraint 
all_distinct. 
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4-1 Binary constraints 

We consider how to implement propagators for the binary constraint A x X = 
B x Y + C, where X and Y are domain variables, A and B are positive integers, 
and C is an integer of any kind. Similar propagators can be implemented for other 
types of binary constraints. 

4-1.1 Forward checking 

Recall that forward checking enforces node consistency. The following shows a prop- 
agator that performs forward checking for the binary constraint. 

<aX=bY+c' (A,X,B,Y,C) => 

'aX=bY+c_forward' (A,X,B,Y,C) . 

'aX=bY+c_forward' (A,X,B, Y,C) , var (X) ,var(Y) ,{ins(X) ,ins(Y)> => true. 
'aX=bY+c_forward' (A,X,B, Y,C) , var (X) => 

T is B*Y+C, X is T//A, A*X=:=T. 
'aX=bY+c_forward' (A,X,B,Y,C) => 

T is A*X-C, Y is T//B, B*Y=:=T. 

The operation opl//op2, which is equivalent to truncate (opl/op2), gives the 
integer quotent of the division. When both X and Y are variables, the propagator 
is suspended. When either variable is instantiated, the propagator computes the 
value for the other variable. 

4-1-2 Interval consistency 

The following propagator, which extends the forward-checking propagator, main- 
tains interval consistency for the constraint. 

c aX=bY+c' (A,X,B,Y,C) => 

'aX=bY+c_reduce_domain' (A,X,B,Y,C) , 
'aX=bY+c_f orward' (A,X,B,Y,C) , 
'aX=bY+c_interval' (A,X,B,Y,C) . 

The subgoal 1 aX=bY+cjreduce_domain' (A,X,B,Y,C) preprocess the constraint 
to make it interval-consistent when the constraint is generated. 

1 aX=bY+c_reduce_domain '(A,X,B,Y,C) => 

'aX in bY+c_reduce_domain' (A,X,B,Y,C) , 
MC is -C, 

'aX in bY+c_reduce_domain' (B,Y,A,X,MC) . 

'aX in bY+c_reduce_domain' (A,X,B,Y,C) => 
L is (B*min(Y)+C) /> A, 
U is (B*max(Y)+C) /< A, 
X in L. .U. 



Programming Finite-Domain Constraint Propagators in Action Rules 13 



The operation opl /> op2 returns the lowest integer that is greater than or equal 
to the quotient of opl by op2 and the operation opl /< op2 returns the greatest 
integer that is less than or equal to the quotient. It can be proved easily that no 
value outside the range L. .U satisfies the constraint. 

The subgoal ' aX=bY+c_interval ' (A,X,B,Y,C) maintains interval consistency 
for the constraint. 

'aX=bY+c_interval> (A,X,B,Y,C) => 

'aX in bY+c_interval' (A,X,B,Y,C) , 
MC is -C, 

'aX in bY+c_interval' (B,Y,A,X,MC) 

'aX in bY+c_interval' (A,X,B,Y,C) , 
var(X) ,var(Y) , 
{generated , bound (Y) } 
=> 

'aX in bY+c_reduce_domain' (A,X,B,Y,C) . 
'aX in bY+c_interval' (A,X,B,Y,C) => true. 

Notice that the action ' aX=bY+cjreduce_domain' (A,X,B,Y,C) is executed only 
when both variables are free. If either one turns to be instantiated, then the forward- 
checking rule takes care of that situation. 

4-1.3 Arc consistency 

The following propagator, which extends the one shown above, maintains arc con- 
sistency for the constraint. 

<aX=bY+c' (A,X,B,Y,C) => 

'aX=bY+c_reduce_domain' (A,X,B,Y,C) , 
<aX=bY+c_f orward' (A,X,B,Y,C) , 
'aX=bY+c_interval' (A,X,B,Y,C) , 
'aX=bY+c_arc' (A,X,B,Y,C) . 

' aX=bY+c_ar c ' ( A , X , B , Y , C) => 

'aX in bY+c_arc' (A,X,B,Y,C) , °/ reduce X when Y changes 
MC is -C, 

'aX in bY+c_arc' (B,Y,A,X,MC) .'L reduce Y when X changes 

'aX in bY+c_arc' (A,X,B,Y,C) ,var(X) ,var(Y) ,{dom(Y,Ey)> => 
T is B*Ey+C, 
Ex is T//A, 

(A*Ex=:=T -> exclude (X, Ex) ; true) . 
'aX in bY+c_arc' (A,X,B,Y,C) => true. 

Whenever an element Ey is excluded from the domain of Y, the propagator ' aX in 
bY+c_arc' (A,X,B,Y,C) is activated. If both X and Y are variables, the propagator 



7o reduce X when Y changes 
. 7. reduce Y when X changes 
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excludes Ex, the counterpart of Ey, from the domain of X. Again, if either X or Y 
becomes an integer, the propagator does nothing. The forward checking rule takes 
care of that situation. 



4-2 Non-binary Constraints 

In indexical-based CLP(FD) systems, constraints are split into indexicals that con- 
tain no more than three variables. This algorithm has several advantages. Firstly, 
it generates linear-size code. Secondly, indexicals can be implemented in a low-level 
language to achieve better performance. Thirdly, information propagation can be 
restricted to only those constraints for which the domains have the possibility to 
be reduced ( |Codognet and Diaz 1996). For example, consider the two ternary con- 
straints Tl = XI + X2 and Tl + A3 + A4 = 0. If Tl = XI + X2 is activated 
by an update of XI, as long as the shared variable Tl does not change the other 
constraint needs not be activated. The disadvantages of this algorithm are that new 
domain variables have to be introduced and the granularity of constraints becomes 
smaller and thus context switching becomes more costly. In B-Prolog, each domain 
variable takes at least 10 words, letting alone the space for the constraints and 
data structures for the elements. The space overhead cannot be neglected when the 
number of variables is large. 

In comparison with indexicals, the high descriptive power of AR opens new ways 
to compiling non-binary constraints. We present two algorithms. One is called unite, 
which adopts one propagator for each constraint to maintain interval consistency. 
The other one, called hybrid, maintains interval consistency when the constraint 
contains more than two variables and maintains arc consistency when the constraint 
turns into binary. 



4-2.1 Unite: use one propagator for each constraint 

Let A\ X X\ + . . . + A n x X n + C = be an n-ary constraint where each Ai 
(i = 1, . . . , n) is a non-zero integer and each Xi is a domain variable or an integer. 
The propagator for the constraint takes the following form: 

'A1*X1+. . . +An*Xn+C=0 ' (C,A1,A2, . . . ,An,Xl,X2, . . ,Xn) , 

{generated, ins (XI) ,bound(Xl) , . . . ,ins(Xn) ,bound(Xn)} 

=> 

... °/« reduce the domains of XI, . . . ,Xn. 

In the action, attempts are made to reduce the lower and upper bounds of the 
domain of every variable. 

To facilitate the generation of the code for reducing domains, the compiler splits 
the expression A\ x X\ + . . . + A n x X n + C into the following list of sub-expressions 
each of which contains at most three variables: 
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T Q = C, 

Ti = T + A 1 *X 1 , 
T 2 = T 1 + A 2 *X 2 , 

T n = T n —i + A n * X n 
The generated reducer first computes the lower and upper bounds of the temporary 
variables 9 by propagating information forward from T to T n . The lower and upper 
bounds of T are computed from those of Tj_i and Ai x Xi (i = 1, . . . , n). After 
that, the reducer propagates information backward from T n to Ti. For each tuple 
Ti = Tj_i + ^ x the new bounds of Tj_i and are computed from those of 
T- 

For example, the following shows the propagator generated for the constraint 
X1+X2+X3+C = 0. 

' Xl+X2+X3+C=0 ' (C , XI , X2 , X3) 

{generated, ins (XI) ,bound(Xl) ,ins(X2) ,bound(X2) , 
ins(X3) , bound (X3)> 

=> 

<X1+X2+X3+C=0_reducer' (C,X1,X2,X3) . 



'Xl+X2+X3+C=0_reducer> (C,X1,X2,X3) => 

Ltl is C+min(Xl), Utl is C+max(Xl) , 7. TI = C+Xl 

Lt2 is Ltl+min(X2) , Ut2 is Utl+max(X2), 7. T2 = T1+X2 
Lt3 is Lt2+min(X3) , Ut3 is Ut2+max(X3) , 7. T3 = T2+X3 
Lt3 =< 0, Ut3 >= 0, 7, T3 = 

% 

NewLx3 is 0-Ut2, NewUx3 is 0-Lt2, 7. T3 = T2+X3 

X3 in NewLx3 . . NewUx3 , 

NewLt2 is 0-max(X3) , NewUt2 is 0-min(X3) , 
7. 

NewLx2 is NewLt2-Utl, NewUx2 is NewUt2-Ltl,% T2 = T1+X2 
X2 in NewLx2 . . NewUx2 , 

NewLtl is NewLt2-max(X2) , NewUtl is NewUt2-min(X2) , 
7. 

NewLxl is NewLtl-C, NewUxl is NewUtl-C, 7. TI = C+Xl 
XI in NewLxl. .NewUxl, 

NewLtl-max(Xl) =< C, NewUtl-min(Xl) >= C. 

The advantage of this algorithm is that only one propagator is used for each con- 
straint whose code size is linear in the number of variables in the constraint. The 
weakness of this algorithm is that the reducer is not fast. Whenever a variable is in- 
stantiated or a variable's bound is updated, the reducer tries to reduce the domains 
of all the variables including the seed variable that triggers the propagator. 



Temporary variables arc plain variables, not domain variables. Therefore, this compilation 
scheme is different from compiling constraints into indexicals. 
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4-2.2 Hybrid: combining interval and arc consistency algorithms 

For a non-binary constraint, it is too expensive to maintain arc consistency. One 
practical strategy is to maintain interval consistency while there are multiple vari- 
ables in the constraint and to maintain arc consistency when the constraint turns 
into binary. The following shows the propagator for the linear non-binary constraint 
AlxXli ... + AnxXn + C=0. 

'A1*X1+. . .+An*Xn+C=0' (C,A1,A2 An,Xl,X2, . . ,Xn) , 

n_vars_gt ( [XI , . . . , Xn] , 2) , 

{generated, ins (XI) ,bound(Xl) , . . . ,ins(Xn) ,bound(Xn)} 

=> 

... 7 reduce domains of XI Xn. 

'A1*X1+. . .+An*Xn+C=0' (C,A1,A2, . . . ,An,Xl,X2, . . ,Xn) => 

nary_to_binary([C,Al,X2,A2,X2, . . . ,An,Xn] ,NewC,Bl ,B2, Yl ,Y2) , 
call_binary_constraint_propagator (NewC , Bl , Yl , B2 , Y2) . 

The propagator is activated whenever any variable is instantiated or its bound is 
updated. When n_vars_gt( [XI, . . . ,Xn] ,2) succeeds, i.e. when there are multiple 
variables in the constraint, the domains are reduced to make the constraint interval- 
consistent. When the constraint becomes binary, the condition n_vars_gt fails and 
the second rule is tried. The subgoal nary_to_binary transforms the constraint into 
the binary constraint Bl xYl + B2 x Y2 + NewC=0, and the next subgoal invokes 
an appropriate propagator for the binary constraint. 10 

4-3 Propagators for all_distinct 

The constraint all_distinct (L) holds if the elements in L are pairwise different. 
One naive implementation method for this constraint is to generate binary dise- 
quality constraints between all pairs of variables in L. This implementation has 
two problems: First, the space required to store the constraints is quadratic in the 
number of variables in L; Second, splitting the constraint into fine-grained ones may 
lose possible propagation opportunities ( |Regin 1994| |Puget 1998D . This subsection 
presents two propagators for the constraint. The propagation algorithms are not 
new. The goal of this subsection is to illustrate the expressive power of AR. 

4-3.1 A linear-space propagator 
To solve the space problem, we define all_distinct in the following way: 
all_distinct (L) => all_distinct (L, [] ) . 

all_distinct ( [] ,Lef t) => true. 

10 In the implementation in B-Prolog, the two built-ins n_vars_gt and nary_to_binary do not take 
the constraint as an argument but instead access the constraint in the parent subgoal. In this 
way, no copy of the constraint needs to be made. 
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all_distinct( [X| Right] , Left) => 
outof (X, Left, Right) , 
all_distinct (Right , [X I Left] ) . 

outof (X, Left, Right) , var(X), {ins(X)} => true. 

outof (X, Left, Right) => exclude_list (X,Lef t) ,exclude_list(X, Right) . 
exclude_list(X, [] ) . 

exclude_list(X, [Y|Ys]) :- exclude(Y,X) ,exclude_list(X,Ys) . 

For each variable X, let Left be the list of variables to the left of X and Right be the 
list of variables to the right of X in L. The subgoal outof (X,Lef t .Right) holds if X 
appears in neither Left nor Right. Instead of generating disequality constraints be- 
tween X and all the variables in Left and Right, the subgoal outof (X , Left , Right) 
suspends until X is instantiated. After X becomes non- variable, exclude_list (X,Lef t) 
and exclude_list(X, Right) exclude X from the domains of the variables in Left 
and Right, respectively. 

There is one propagator outof (X, Left , Right) for each element X in the list, 
which takes constant space. Therefore, all_distinct (L) takes linear space in the 
size of L. Notice that the two lists Left and Right are not merged into one bigger 
list; otherwise, the constraint would take quadratic space. 

4-3.2 Weak arc consistency 

In terms of pruning ability, the linear-space propagator is the same as the naive 
one that splits a constraint of all_distinct into binary disequality constraints. In 
this subsection, we present a propagator that has stronger prunning power than 
the naive propagator. 

Given any set of values D of size n, the constraint all_distinct (i) is said to be 
weak arc consistent if there are at most n variables in L whose domains are subsets 
of D. For each variable X in L, let L — {X} be the list of variables in L but X, n be 
the size of the domain of X, and m be the number of variables in L — {X} whose 
domains are subsets of that of X. If m + 1 > n, then the constraint is unsatisfiable 
since it is impossible to assign n values to more than n variables such that each 
variable gets a different value. If m + 1 = n, then for each value v in X's domain, 
we can safely exclude v from the domains of all the variables whose domains arc 
not subsets of that of X. 

Consider the following query, 

X in {1,2}, Y in {1,2}, Z in {1,2}, all_distinct( [X,Y,Z] ) . 

the weak arc-consistency propagator detects the inconsistency of the constraint 
without labeling any variables. For the following query, 

X in {1,2}, Y in {1,2}, Z in {1,2,3}, all_distinct ( [X, Y,Z] ) . 

the algorithm assigns 3 to Z without labeling any variables. 
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The weak arc-consistency propagator is not as powerful as the algorithm proposed 
by Regin (Reg in 1994| ) in terms of pruning ability but is much easier to implement. 
To incorporate weak arc-consistency checking into the linear-space propagator, we 
only need to redefine outof (X, Left .Right) as follows: 

outof (X, Left .Right) , var(X), {generated, ins(X) ,bound(X) ,dom(X)} => 

outof _reducer (X, Left , Right) . 
outof (X,Lef t , Right) => exclude_list (X,Lef t) , exclude_list (X, Right) . 

where outof _reducer(X, Left , Right) first counts the variables in Left and Right 

whose domains are subsets of the domain of X and then decides what action to take 
depending on the count and the size of the domain of X. 

The key operation is to decide whether a domain is a subset of another domain. 
In the worst case, the two domains have to be scanned. There are several facts that 
can be used to avoid scanning domain elements. A domain Dl cannot be a subset 
of another domain D2 if Dl has a larger size or has a larger interval. Also if two 
domains are intervals without holes, then scanning the elements is unnecessary. 
Another fact that can be used in the detection is that if the event is dom(X,E) 
meaning that E has been excluded from X's domain, then another domain Y cannot 
be a subset of X if E is included in Y. To take advantage of this fact, the propagator 
can be rewritten into the following: 

all_distinct (L) => all_distinct (L, [] ) . 

all_distinct ( [] ,Lef t) => true. 
all_distinct( [X| Right] .Left) => 

outof (X , Lef t , Right ) , 

outof _dom(X, Left , Right) , 

all_distinct (Right , [X I Left] ) . 

outof (X, Lef t .Right) , var(X), {generated, ins (X) ,bound(X) } => 

outof _reducer(X, Lef t , Right) . 
outof (X,Lef t .Right) => exclude_list (X,Lef t) , exclude_list (X, Right) . 

outof _dom(X, Left, Right) ,var(X) , {dom(X,E)} => 

outof _reducer(X,E, Left, Right) . 
outof _dom(X, Lef t , Right) => true. 

The subgoal outof _reducer(X,E, Lef t .Right) takes E into account when detecting 
whether a domain is a subset of that of X. 

5 Performance Evaluation 

B-Prolog has been extended to accommodate AR and the finite-domain constraint 
solver described in this paper has been developed in AR. In this section, we evaluate 
the performance of the finite-domain constraint solver. 
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Table 1. Comparison of CPU times. 
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EP 
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1 
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1 
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1 
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1 
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1 
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3.57 


1 





92 
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1 


1 


00 


1 


67 


7.03 


4.67 


1 


1 


60 


color 


1 


1 


14 


1 


00 
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3.01 


1 


1 


25 


eqlO 


1 





98 


3 


77 


4.77 


4.92 


1 


3 


78 
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1 


1 


06 


2 


00 


4.23 


3.34 


1 


1 


96 


magic3 


1 


1 


38 


1 


98 


8.44 


4.41 


1 


1 


52 


magic4 


1 


1 


18 


1 


96 


8.45 


6.27 


1 


1 


57 


Olympic 


1 





75 


2 


25 


11.25 


4.75 


1 


1 


43 


queens(25) 


1 


1 


01 





43 


4.24 


5.03 


1 





43 


sendmoney 


1 


1 


09 


3 


65 


6.74 


7.78 


1 


2 


62 


sudoku81 


1 


1 


00 


2 


28 


6.67 


6.18 


1 


1 


36 


zebra 


1 


1 


13 


2 


33 


7.86 


9.61 


1 


1 


77 



<arithmetic mean> 1 1.04 1.96 6.68 5.13 1 1.61 

<geometric mean> 1 1.03 1.72 6.36 4.84 1 1.42 



We compared the performance of B-Prolog version 6.7 (BP) 11 with three other 
CLP(FD) systems: ECLiPSe 5.8 #77 (EP), GNU-Prolog version 1.2.16 (GP), and 
SICStus 3.12.0 (SP). There are two solvers available in BP: one is called BP-AC 
which adopts the hybrid algorithm presented in this paper for equality constraints 
and the other called BP-IC which maintains only interval consistency for equality 
constraints. BP- AC is the default solver. 12 . The linear-space propagator is used for 
all_distinct in both solvers. 

Table ^ shows the CPU times taken by the four solvers to run a set of bench- 
marks, 13 assuming the time taken by BP-IC be 1. Most of the benchmarks have been 
widely used by other authors to compare CLP(FD) systems IjCarlsson et al. 19971 
ICodognet and Diaz 19961 |Hentenryck 19891 . Three new programs were added by 
the author into the set: color is a program that colors a map with 110 regions; 
Olympic is a puzzle taken from a Mathematics Olympic game for elementary stu- 
dents; and sudoku81 is a program for solving a puzzle. The left-to-right labeling 
strategy is used to instantiate variables in all the benchmarks. The CPU times were 
measured on a 1.7GHz CPU running Windows XP. Each program was run at least 
10 times and the average was taken. For some programs, execution was repeated 
up to 1000 times to obtain a stable average. Garbage collection was disabled. GP 
has a native code compiler for Linux. The comparison with GP was also conducted 
on Linux. 



Available from www.probp.com. 

To switch to BP-IC, set the Prolog flag constr.consistency to int 
Available from probp.com/bench.tar.gz. 
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Table 2. Comparison of numbers of backtracks. 



Program BP-AC BP-IC GP 



alpha 


4605 


8440 


8440 


bridge 











cars 


53 


53 


34 


color 
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560 


560 


eqlO 


49 


49 


49 


eq20 


49 


49 


49 


magic3 


2 


2 


2 


magic4 


18 


18 


18 


Olympic 


36 


50 


50 


queens (25) 


7255 


7255 


7255 


sendmoney 


2 


2 


2 


sudoku81 











zebra 


2 


2 


2 



The BP solvers compare favorably with GP and are significantly faster than EP 
and SP. EP is the slowest among the compared solvers, probably because of the 
overhead of supporting priority-based scheduling ( Wa llace et al. 2004|l . BP outper- 
forms GP remarkably for programs that contain non-binary equality constraints, 
such as eqlO, eq20, and sendmoney. This result reveals that the disadvantages of 
splitting n-ary constraints into indexicals outweigh the advantages. On the other 
hand, GP is more than twice as fast as BP for queens (25). The high speed for 
queens may be attributed to an optimization technique adopted in GP that com- 
bines indexicals. If the propagators for the disequality constraints were combined 
for the program in BP, 14 the speed would be doubled. 

Comparing BP- AC and BP-IC reveals that the hybrid algorithm is effective for 
alpha and Olympic only. The BP-AC solver is adopted as the default one since for 
some programs, such as the queens program given in (Pu get and Leconte 1995| l, 15 
BP- AC is exponentially faster than BP-IC. BP- AC is slightly slower than BP-IC for 
some of the programs. In general, this happens for programs for which the efforts 
to reduce domains do not pay off. 

Table |U gives the numbers of backtracks performed by the three solvers. BP-AC 
makes the same number of backtracks as BP-IC except for alpha and Olympic, and 
GP makes fewer backtracks than BP-AC for cars. Basically, the three solvers explore 
the same search trees for most of the programs. Therefore, the comparison results 
shown in Table reflect the real performance of the solvers. 

GP and BP are quite different. In GP constraints are compiled into indexicals 



14 For the three disequality constraints X Y, X Y + N, and X Y — N , we can use one 
propagator rather than three to handle the ins(X) and ins(Y) events. 

15 This program is not included in the benchmark set since it requires support of negative integer 
domains and thus cannot run on GP. 
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defined in C while in BP constraints are compiled into propagators defined in action 
rules. Although the GP Prolog engine may not have much impact on the perfor- 
mance of constraint programs, the BP engine does have a great impact since all the 
propagators are defined in action rules. One evidence for this observation is that the 
BP constraint solver becomes 20-30 percent faster after the main switch statement 
in the emulator is changed to a jump table. A further speed-up is expected if a 
native code compiler is employed. 

There are other factors that affect the performance of a solver, such as domain 
representation, interaction with other solvers, and garbage collection! Wallac e et al. 20 04). 
GP supports only finite domains of positive integers, while BP supports not only 
finite integer domains but also trees and finite domains of ground terms and sets 
l|Zhou 2002(1 . In BP, integer domains are represented as described in Subsection l2.3l 
BP adopts a sound and complete arithmetic that guarantees that solutions found 
are correct and no solution is lost. When excluding an inner element from a large 
interval domain, the system generates a disequality constraint rather than brutally 
changes the interval into a bit vector as is done in GP. 16 BP has a garbage col- 
lector that collects garbage on the heap and the control stack, but GP does not 
support garbage collection yet. Garbage collection may suppress some optimization 
techniques. 

6 Related Work 

CLP(FD) systems have undergone an evolution process, from closed to open and 
from low level to high level. Several constructs have been proposed to facilitate the 
implementation of constraint propagators. Examples include attributed variables 
(Holzbaur 1992), indexicals ( |Codognet and Diaz 19961 ), extended indexicals called 
projection constraints (ISidebottom and Havens 1996|l . delay clauses iMeier 19941 
IZhou 1 998)1 . and constraint handling rules l|Fruhwirt h 199811 . An action rule is an 
extension of a delay clause that allows for the descriptions of not only delay condi- 
tions on subgoals but also activating events and actions. This section compares AR 
with these constructs introduced into Constraint Logic Programming. Constructs 
introduced into other languages such as ILOG ( |Puget and Leconte 1995) ) and Oz 
(Schultc 2002) are not compared. 

AR is more powerful and flexible than indexicals. We have described in this 
paper several propagation algorithms in AR, some of which cannot be encoded 
in indexicals (e.g., the hybrid algorithm for n-ary constraints) and some of which 
cannot be implemented as efficiently (e.g., arc consistency for binary constraints). 
Consider the following indexical taken from ( Carlsso n et al. 19 97). 

X in dom(Y) + C 

which maintains arc consistency for the constraint X = Y + C w.r.t. X. Whenever 
an element y is excluded from the domain of Y, the indexical is activated. Because 



Generating a disequality constraint is less efficient than changing an interval into a bit vector 
since the disequality constraint needs to be checked each time the variable is instantiated. 
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the indexical does not know what the excluded element is, it has to go through the 
domain elements of Y in the worst case to locate a possible no-good value in the 
domain of X. In contrast, in the propagator implemented in AR, the propagator 
knows exactly what element is excluded from the domain of Y and thus can compute 
the counterpart in the domain of X in constant time. 

Compiling constraints into indexicals enables the use of more specialized propaga- 
tors and restricts propagation to within a small number of constraints fCodognet and Diaz 1996] 
Nevertheless, this approach has to introduce new temporary variables and lower 
the granularity of propagators. Our experiment reveals that B-Prolog outperforms 
GNU-Prolog for almost all the benchmarks that contain non-binary constraints. 
This result reveals that the disadvantages of splitting constraints outweigh the ad- 
vantages. Similar observations have been made independently in f Oar lsson et al. 19971 
|Harvey and Stuckey 2003|IZhou 1998| . In flHarvey and Stuckey 2003| l, the same two- 
phase algorithm is used to reduce domains of linear constraints. 

Attributed variables IHolzba ur_1992|| are variables with attached attributes each 
of which has a list of handlers. Touching an attribute triggers the corresponding 
list of handlers. In order to make context switching swift for handlers, systems such 
as ECLiPSe treats handlers as demons rather than as normal subgoals. A demon 
is different from a normal subgoal in that it does not disappear after execution but 
instead waits for another activation. In this sense, agents in our system are similar 
to demons. Nevertheless, an agent can be activated by different kinds of events and 
an agent may take different actions depending on the conditions. An agent can 
be defined by multiple action rules and the rules are compiled into a tree by the 
compiler such that shared tests are combined and conditions that failed once need 
not be tested again. SICStus l|Carlsson et al. 1997(1 provides interfaces for imple- 
menting propagators and also some sort of delay construct similar to attributed 
variables that triggers propagators when events are posted. 

AR is an extension of our early delay construct proposed in IjZhou 1998fl that 
allows for the event dom(X,E) and user-defined events. The support of the event 
dom(X,E) is essential for implementing arc consistency algorithms and also prop- 
agators for set constraints (Zhou 2002). Our delay construct is an extension of 
Meier's delay clause construct IMeier 1994|l that allows for not only delay condi- 
tions but also events and actions. In Meier's delay clause, events are implicitly 
extracted from delay conditions and a delayed subgoal never takes actions as long 
as the delay condition is satisfied. In retrospect, all these constructs were inspired 
by early work by Colmerauer and Naish IjColmerauer 19841 INaish '1 985 ) . 

Other rule-based languages have been designed for implementing constraint prop- 
agators. CHR resembles a production system. In CHR, the left-hand side of a rule 
specifies a pattern of constraints in the constraint store and the right-hand side 
specifies new constraints to replace those on the left-hand side or to be added into 
the store. It should be possible to implement in CHR all the propagation algorithms 
described in this paper provided certain built-ins are added. If events are treated as 
constraints, then an action rule can be translated into a CHR rule. Treating events 
as constraints, however, can hardly achieve the same performance. Events are re- 
moved automatically after all the agents that are waiting for them are activated. 
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In CHR, there must be rules to remove the events explicitly. The left-hand sides 
of CHR rules can have multiple constraint patterns. Therefore, it is impossible in 
general to translate a CHR rule into action rules straightforwardly. It is not clear 
whether or not it is possible to simulate CHR rules in action rules and how if the 
answer is yes. It would be an interesting direction to explore in the future. 

7 Conclusion 

There is a need for an implementation language for constraint propagators that 
is expressive enough and can be implemented efficiently This paper has presented 
such a language called AR. The expressiveness of the language is illustrated though 
several examples that cannot be implemented in indexicals: the propagator for 
maintaining arc consistency of binary equality constraints; a weak arc-consistency 
propagator for the all_distinct constraint; and a hybrid algorithm for non-binary 
equality constraints that combines interval and arc consistency ones. The efficiency 
is evaluated through benchmarking. For a set of widely used benchmarks, our solver 
implemented in B-Prolog is significantly faster than that of GNU-Prolog, one of the 
fastest finite-domain constraint solvers available now. 

The results are encouraging and promising since our solver is implemented in 
a high-level language and B-Prolog is an emulator-based system which provides 
more facilities than GNU-Prolog such as garbage collection and constraint solving 
over other domains. The high-performance of our solver stems from the follow- 
ing facts. Firstly, only one propagator is generated for each non-binary equality 
constraint that maintains interval consistency. Our solver performs especially well 
for the benchmarks that contain non-binary equality constraints. This reveals that 
compiling non-binary equality constraints into indexicals has more cons than pros. 
Secondly, the hybrid algorithm adopted in our solver is a good compromise be- 
tween the need to achieve high-level consistency to cut search spaces and the need 
to reduce the cost. The cost of achieving arc consistency for binary constraints 
is relatively small, but the effect can be very big for certain programs. Thirdly, 
our solver employs optimization techniques that reduce redundant activations of 
propagators. 

Our solver can be improved further in the following aspects: (1) develop new 
optimization techniques for further avoiding redundant activations of propagators; 
and (2) implement consistency algorithms beyond interval and arc consistency such 
as path consistency and Regin's algorithm for all_distinct. 
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